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ABSTRACT 


Pump cavitation compliance represents the hydraulic compressi- 
bility of the bubbles formed from the complex cavitation process at pump 
inlets. Such ,T bubble n compliance is a significant factor in establishing the 
resonant frequencies for the dynamics of rocket engines and, as such, must 
be known for the pogo stability analysis of a launch vehicle. 

Cavitation compliance has been determined experimentally from 
pulsing tests on a number of rocket turbopumps. The primary test data 
used for this study are those for the Rocketdyne H-l, F-l, and J-Z oxidizer 
and fuel pumps employed on Saturn vehicles. The study shows that these 
data can be correlated by a particular form of nondimensionalization, the 
key feature of which is to divide the operating cavitation number or suction 
specific speed by its value at head breakdown. A previous study (Ref. 3) 
concluded that such a correlation was not possible since breakdown infor- 
mation was not employed. An expression is obtained for a best-fit curve 
for these data. Another set of test data for the Aerojet LR87 and 91 pumps 
can be correlated by a somewhat different nondimens ional pump performance 
parameter, specifically by relating the cavitation number to its position 
between the head breakdown point and the point of zero slope of the head 
coefficient versus cavitation number. 

A study of mathematical models for the cavitation compliance 
correlation shows a degree of success. Involved are an assumption that the 
bubble volume at any pump operating point is proportional to a power of a 
nondimensional performance parameter involving a relationship to head 
breakdown. 

Recommendations are given for the estimation of the cavitation 
compliance for new designs in the Rocketdyne family of pumps. 
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SYMBOLS" 


A 



c 
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Cu 

c 

D 

d 


g 

h 

K 

K I 

L 

M 

N 



n 


NPSH 


P 

s 


V 


2 2 
area, m , in. 

2 2 

cavitation compliance, m , in. 

backflow absolute velocity, m/sec, ft/sec 

axial component of the absolute velocity, m/sec, ft/sec 

whirl component of absolute velocity, m/sec, ft/sec 

length of vapor cavity, m, in. 

impeller diameter, m, in. 

blade spacing, 2 irr/z, m, in. 

2 2 

gravitational constant 9*807 m/sec , 32.2 ft/sec 
cavity height at cavity closure, m, in. 
cavitation number 

cavitation number at zero slope of 4 1 versus K curve 

ratio of backflow to through flow 

2 

mass, Kg, lb -sec /ft 

pump rotational speed, revolution per sec, revolution per min (rpm) 

, 3 /4 

pump specific speed, dimensionless, rpm v /gpm/(ft) 
empirical exponent 

Net Positive Suction Head, Joules /Kg, ft 

static pressure, Newton/m^, psia 

2 

vapor pressure. Newton/m , psia 


In describing each symbol the unit immediately following the description is 
in SI units followed by the conventional engineering units. 
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SYMBOLS (Continued) 


3 3 

Qq volumetric flow, m /sec, ft /sec 

R gas constant. Joules /Kg °Kelvin, ft-lb/lb°R 

r radius, m, in. 

r. lower radius, backflow region, m, in. 

S suction specific speed, dimensionless, rpm v /gpm/(ft)^ ^ 

T temperature, °Kelvin, °Rankine 

t stay time of the bubble, sec 

U tangential velocity, m/sec, ft/sec 

3 3 

V volume, m , in. 

W relative velocity, m/sec, ft/sec 

mass fraction of gas in liquid at the pump inlet 

z number of blades, inducer 

a inducer inlet incidence angle, a = (3 - (3^,, Radian, deg 

P inducer inlet blade angle with respect to tangent, Radian, deg 

P-., inducer inlet relative velocity angle, Radian, deg 

r 

y the angle measured between the tangent and the tip circular arc 

at the impeller discharge and the axis of rotation, Radian, deg 

A diffusion factor 

Ah pump head rise. Joule /Kg, ft 

A^ head loss coefficient 

e impeller diameter ratio, r ^/ r 

q pump efficiency 
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SYMBOLS (Continued) 


p. slip factor 

3 3 

p density. Kg/m , lb /ft 

0 backflow velocity absolute angle, Radian, deg 

0 flow coefficient 

t solidity 

2 

4j head coefficient, Ah/(U^/g) 

^ total head coefficient, 47^ 

4j actual Euler head coefficient 

Th 

Jj ideal Euler head coefficient 

oo 



b bubble 

c quantities evaluated on free streamline 

Ext External 


G gas 

h inducer inlet hub 


i inducer inlet eye 

j any station along the radius 

L liquid 

1 inducer inlet, upstream 

2 inside inducer at wake collapse 


-xi- 



SYMBOLS (Concluded) 


m mean 

s static 

T total 

t impeller tip 

Superscript 


at cavitation breakdown 



1. INTRODUCTION 


1. 1 BACKGROUND 

Longitudinal structural oscillations and corresponding propellant 
oscillations at frequencies in the range of 5 to 60 Hz, referred to as pogo, 
have been observed to erupt spontaneously on a number of liquid propellant 
launch vehicles during powered flight (Ref. 1). Pogo results from an unstable 
dynamic interaction of the vehicle (structure and tanks) with the propulsion 
system (engines and feedlines). The propulsion system displays resonant 
characteristics as a result of hydraulic resonances of the feedlines. The 
relationship of the frequencies of feedline resonances to the natural fre- 
quencies of the vehicle structure is a major factor for system stability. A 
feedline resonant frequency can be significantly affected by (1) pump cavita- 
tion compliance, (Z) the structural or hydraulic compliance of devices such as 
accumulators or pressure -volume compensators, and (3 ) the distributed 
compliance along the feedline due to propellant compressibility and pipe-wall 
radial elasticity. Cavitation compliance is defined as 

r dV ” 

C b ' dP P L g 

S 

L.e., the product of the volumetric compliance (negative rate of volume 

Z 

change per unit pressure) and propellant weight density; its units are meters , 

Z 

and in engineering units, inches . Pump cavitation compliance represents 
the hydraulic compressibility associated with bubbles formed from the com- 
plex cavitation process at pump inlets. These bubbles have been observed in 
flow- visualization studies, and this bubble compliance must be known for the 
pogo stability analysis of a launch vehicle. 

Pump cavitation compliance could be determined directly by meas- 
uring oscillatory flow both upstream and downstream of the cavitation region 
in conjunction with oscillatory pump-inlet pressure measurement. This has 
generally not been possible because of the unavailability of dynamic flowmeters. 
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As a result the cavitation compliances must be inferred via a less direct 
means* Often the inference is based on the pump-inlet compliance required 
to yield an observed resonant frequency based on a mathematical model of 
the particular test or flight configuration. Obviously, the result can reflect 
errors resulting from any misrepresentation in the system. All the cavita- 
tion compliance data used in this study is subject to uncertainty resulting 
from such error. Some of these uncertainties are noted and discussed in 
Ref. 3 For purposes of this study, only the best estimates of the cavita- 
tion compliances are employed. 

An earlier empirical study of the cavitation compliance of Aerojet 
pumps used on Titan vehicles was conducted by Ghahremani (Ref. 2). That 
study concluded that the compliance is primarily the result of backflow 
(recirculation to the pump inlet of high-pressure flow past the tip of the 
unshrouded impeller) and not the result of cavitation actually attached to the 
inducer blades. The result of that study was a partial agreement between 
an empirical model and the experimental data. 

A more recent study was conducted by Martin Marietta Corporation 
(Ref. 3) using the same Saturn vehicle pumps and Titan pumps that are con- 
sidered in this study. The main emphasis of that study was to develop an 
analytical model using a potential flow solution to determine cavitation com- 
pliance by calculating vapor volume between the pump blades caused by blade 
cavitation. Gaseous cavitation, and cavitation resulting from flow recircu- 
lated to the pump inlet were not treated. A computer program was developed 
and applied to the following oxidizer pumps: J-2, F-l, and H-l Saturn pumps 

and the LR87 Titan pump. The analytical method predicted much less cavi- 
tation compliance than the test data indicates. The results were quite dis- 
appointing, with no cavitation (or practically none) predicted for the F-l, 

Not mentioned in Ref. 3 is the possible importance of pump dynamic gain 
in deriving cavitation compliance. In a current study of the Thor/Delta 
Vehicle and its MB-3 engine, we have recently found that pump dynamic gain 
(partial derivative of head rise with respect to inlet pressure, Ref. 1) can 
substantially influence the resonant frequency of a feedline /engine system. 

It would be worthwhile to review the various determinations of cavitation 
compliance to determine how the pump gain was considered. 
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H-l, and LR87 pumps. It was concluded that mechanisms other than blade 
cavitation must predominate. An empirical study of the data-in Ref. 3, 
led to the additional conclusion that a nondimensional combination of pump 
geometry and fluid properties could not be found to correlate the available 
test data. 

1.2 OBJECTIVE 

The objective of this study was to develop an empirical method which 
relates the F-l, J-2, and H-l pump-inlet cavitation compliance to turbo - 
pump design geometry, performance characteristics and fluid properties. 

The empirical model, as developed, could be used for the prediction of the 
cavitation compliance for new designs of this class of pumps in advance of 
experimental determinations. Such prediction should improve the quality 
of preliminary pogo stability analysis of new vehicles such as the Space 
Shuttle . 

1. 3 SCOPE 

This study was primarily concerned with the prediction of cavitation 
compliance for the F-l, J-2, and H-l pumps operating at their normal speed 
under flow conditions for which the compliance data were available. The 
study considered normal blade cavitation, backflow cavitation, and their 
effect on compliance. In addition to the Saturn pumps. Titan pumps were 
also studied as an aid to understanding pump design effects. 
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2. EVALUATION OF DATA 


2. 1 CAVITATION COMPLIANCE TEST DATA 


The cavitation compliance data as a function of cavitation number 
at the inducer eye is shown in Fig. 1 and is based upon Fig. 4-2M of Ref. 3. 
Cavitation number at the inducer eye is a dimensionless parameter defined 
as 


K. 

i 




( 2 - 1 ) 


where 



(2-2) 


Substituting Eq. (2-2) into Eq. (2-1), we obtain 


K. 

i 


NPSH - 



2 



(2-3) 


Assuming no pre whirl at the reducer inlet 

? 2 2 

W = U + C 

i l m 


(2-4) 


Omitted is data for the MB -3 pumps. Our current investigations of pogo on 
Thor/Delta vehicles have shown that the MB-3 fuel pump data is invalid, 
and the oxidizer pump data requires further study. 

To comply with the contractual requirements of this study, all the deriva- 
tions are in SI units unless otherwise specified. 
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Fig. 1. Oxidizer and Fuel Pumps Cavitation Compliance 




and employing the flow coefficient 


0 


1 


C 


m 

U. 


1 


we obtain from Eq. (2-3) the desired form 


K. 

l 


/npsh 


V 


U?/2 



* 


(2-5) 


For use in our studies, modifications of the F-l oxidizer, J-Z oxidizer, and 
H- 1 fuel data are made as shown by the dotted curves. These modifications 
are discussed in the following text. 

2. 1. 1 F-l Oxidizer Pump Compliance Data 

The solid curves in Fig. 2a (from Fig. 4. 9 of Ref. 3) display uncer- 
tainty in the cavitation compliance of the F-l oxidizer pump resulting from 
differences in alternative determinations of the compliance of a pressure- 
volume compensator (PVC). The corresponding PVC compliances are shown 
in Fig. 2b (from Fig. 4. 7 of Ref. 3). The cavitation compliance is deter- 
mined by subtracting the PVC compliance from the overall compliance as 
inferred from an observed resonant frequency during operation of the pump. 
On the one hand, non-flow feedline dynamic tests indicate a PVC compliance 
which decreases rapidly with static pressure. However, static pressure 
testing of the PVC itself indicates a more constant compliance with pressure, 

although compliance is reduced by mechanical nonlinearity at pressures above 
2 

0. 7 MN/m absolute (100 psia); see Figs. 4. 4 and 4. 5 of Ref. 3. It is con- 
cluded in Ref. 3 that the best estimate of PVC compliance lies between the 
variable and contant estimates (see dotted curve in Fig. 2b), In our opinion 

the PVC compliance, shown by the dashed curve in Fig. 2b, is a more 

reasonable estimate and this is the basis for the dashed curve in Fig. 2a, 

and for the dashed curve of the F-l oxidizer pump in Fig. 1. 
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0.4 0.6 0.8 1.0 
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Fig. 2. F-l Oxidizer Pump Cavitation Compliance and PVC 
Compliance Versus Pressure 
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2 . 1 . 2 


J-2 Oxidizer Pump Compliance Data 


Figure 4. 16 of Ref. 3 shows a large discrepancy in the results of 
independent evaluations of the test data from several facilities to derive 
cavitation compliance. Brown Engineering employed a single -compliance 
pump model to analyze data from S-II and S-IVB feed system tests. Rocketdyne 
employed a double -compliance model to analyze their test facility results. 
Reference 3 states that the results are consistent insofar as resonant fre- 
quency is concerned, but produce different cavitation compliance results 
because of the two types of pump models. Reference 3 rates the Rocketdyne 
model as the better one, but states that the true cavitation compliance can lie 
anywhere between the two results. In our conversion of the Brown and 
Rocketdyne compliances versus inlet static pressure, to versus cavitation 
number, we find that the curve in Fig. 4. 21 of Ref. 3 is nearly equal to the 
Brown results rather than being an intermediate result. For our purposes, 
we have elected to use the more reliable Rocketdyne result and this is shown 
by the dashed curve for the J-2 oxidizer pump in Fig. 1. 

2. 1. 3 H-l Fuel Pump Compliance Data 

Reference 3 employs an erroneous value for the inducer -eye dia- 
meter of the H-l fuel pump, using instead of the large inlet diameter the 
smaller diameter at the tip of the inducer discharge. The dotted curve for 
the H-l fuel pump in Fig. 1 is based on the use of the correct diameter to 
establish the cavitation number. 

2.1.4 LR91 Fuel Pump Compliance Data 

Compliance data for the LR91 (Titan Stage II) fuel pump, obtained 
from feedline resonance observed in flight, is limited to a single NPSH. It 
appeared initially that recent static firing tests would provide Stage II fuel 
pump data over a wide range of NPSH and would extend the oxidizer pump 
data over a much wider range of NPSH than available from flight. Oscillo- 
graphic pump-inlet pressure data were obtained from Aerojet and reviewed 
visually in an attempt to identify feedline resonant frequency versus NPSH. 
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Unfortunately, feedline resonant frequencies could not be detected and this 
search for additional Titan data was discontinued. 

2. 2 PUMP DESIGN AND PERFORMANCE DATA 

Two sources of pump performance and design data were used in this 
study. Initially, the work in Ref. 3 was available and then further data were 
obtained from Rocketdyne. The Rocketdyne data are presented in Table I, 
together with similar Titan data which were obtained from Refs. 2 and 4. 

The pump data in Ref. 2 were based on pump tests without pre valves up- 
stream of the pumps, whereas the data shown in Ref. 4 were based on pump 
tests with prevalves. Pre valves are used in the flight vehicle and they have a 
substantial effect on the pump performance because of the high backflow of 
these pumps. 

In addition to the data shown in Table I, the leading-edge blade -angle 
distributions for F-l, J-2, and H-l pumps were received from Rocketdyne 
and are presented in Fig. 3. For the Titan pumps, it is assumed that the 
leading-edge blade -angle distribution is that of a radial element 

r tan (3 = Constant (2-6) 

where r is the radial distance from axis of rotation and (3 the leading-edge 
blade angle. 

2. 2. 1 Cavitation Performance 

In general, the dimensionless cavitation performance of a pump can 
be presented in terms of head coefficient (4 j ) versus cavitation number (K) 
for a fixed flow coefficient (0^) as shown in Fig. 4. At high K values, the 
slope of the curve is zero and as K is reduced to Kj, the slope begins to 
change. As K is reduced further to K*, the head breaks down and the slope 
of the curve becomes infinite. The cavitation performance data were used 
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Table I. Design and Performance Data for Rocket Engine Pumps 



10 


Notes: ill Performance data apply to the test conditions for which cavitation complu 

determined. 

(ill All pumps are centrifugal except for the axial J-2 fuel pump. 




o 

TT 



U! U313WVI 


Inducer Leading-Edge Blade Angle Distribution for 
Rocketdyne Pumps 
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to establish the he ad -breakdown cavitation number (K*) for each pump as a 
function of the pump-inlet flow coefficient (0^) where 


K* 



(2-7) 


0 


1 


= C /U. 
m l 


( 2 - 8 ) 


where NPSH* = Net Positive Suction Head at cavitation head breakdown in 
m(ft). The performance curves received from Rocketdyne^ and Ref. 4 were 
used for the calculation of K*. The data are shown in Figs. 5 and 6 for the 
Rocketdyne and Titan pumps, respectively. Figure 5 shows the results of 
several tests using inducers alone, model pumps, and full-scale pumps. All 
were tested with water, except for the full-scale J-Z fuel pump. All full- 
scale pumps were tested with propellant, except for the H-l fuel pump. 

The NPSH* values used to calculate K* were determined from pump 
performance data such as that presented in Fig. 7, which shows the cavita- 
tion performance of the F-l oxidizer pump. For a value of about 95 percent 
of the non-cavitating head, the slope is practically infinite. Figure 5 shows 
the K* values determined at the point where the slope of head versus NPSH 
is essentially infinite. For the F-l, J-2, and H-l pumps this occurs at 
values between 90 and 95 percent of the non-cavitating head, the precise 
value being difficult to establish because of the extreme sensitivity to data 
accuracy. For uniformity, Fig. 5 is based on 95 percent of non-cavitating 
head for these pumps. Again for uniformity, 90 percent of non-cavitating 
head for the Titan pumps was arbitrarily selected to determine the K * 
shown in Fig. 6, since data were not available at larger head loss for one 
of the pumps. When data were available, the head loss required to reach 
infinite slope for the Titan pumps was about 80 percent. 

*r ' 

3 Personal communication by letter 71RC6543 from P. N. Fuller to 
F. G. Ghahremani, dated 2 September 1971 
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Fig. 5. Cavitation Number at Breakdown for the Rocketdyne Pumps 




LR87 OX 
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16 - 


NET POSITIVE SUCTION HEAD (NPSH), ft 





The difference between the cavitation performance with water and 
the propellant is due to an effect called the "thermodynamic suppression 
head, 11 (TSH). In general, the formation of vapor bubbles requires heat from 
the surrounding liquid. Therefore, the liquid is cooled as bubbles are formed 
and the vapor pressure of the liquid is reduced. This depression of the vapor 
pressure is called TSH. Liquid hydrogen has the highest TSH value (90 to 
130 ft of head). Liquid oxygen has values of the order of 5 to 15 ft. RP-1 
and Aerozine are known to have negligible TSH. There are indications that 
might have some finite TSH value, but no tests have been performed 
for its determination. 

The results shown in Fig. 6 are for the Titan pumps tested with 
water; both the fuel (Aerozine -50) and the oxidizer (N^O^) are assumed to 
have zero TSH. Note that for the F-l fuel pump several tests were performed 
with the inducer alone, with model pumps using water, and with the full- 
scale pump using both water and RP-1; note also that both the model F-l 
oxidizer and fuel pumps tested with water have a higher cavitation number 
at breakdown than do the full-scale pumps. Furthermore, observe that the 
F-l fuel pump tests with water and RP-1 show about the same cavitation per- 
formance . 
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3. EMPIRICAL CORRELATION OF CAVITATION 
COMPLIANCE DATA 


Our attempts to collapse the pump cavitation compliance data by 
use of pump performance and geometrical information were based on the 
compliance data presented in Fig. 1, using the dashed curves for the F-l 

oxidizer, J-2 oxidizer, and H- 1 fuel pumps. These basic data are presented 

2 2 

as cavitation compliance in m (in. ) versus cavitation number K, Two 

unsuccessful attempts were made in Ref. 3 to collapse the curves by non- 
dimens ionalizing the compliance data. In one case, the compliance was 

divided by the number of inducer blades times the square of the inducer -eye 

2 

diameter, D- (see Fig. 4.22 of Ref. 3). In the other case, the was 
^7 7 

replaced by (D^ - D^), where is the inducer hub diameter, and the recip- 
rocal of the suction specific speed was employed in place of the cavitation 
number K (see Fig. 4. 23 of Ref. 3). In our study we have elected to non- 
dimens ionalize the cavitation compliance by dividing by the pump-inlet area. 
Moreover, while K is a dimensionless cavitation parameter, other pump 
performance information is also required to characterize the degree of 
cavitation. For example, the relationship of the pump operating point to the 
head breakdown point must also be important. As our first attempt we take 
the ratio of K/K* to indicate the cavitation status at the operating point. 

Figure 8 shows the dimensionless compliance C^/A^ versus K/K* for all the 
pumps. With the exception of the J-2 fuel pump, the curves for the Rocketdyne 
pumps collapse rather well. The curves for the Titan pumps exhibit a much 
steeper slope than do those for the Rocketdyne pumps and are not collapsed 
by the plotting. 

A second possibility is to plot the data versus the ratio S/S* (ratio 
of suction specific speeds). In turbomachinery engineering practice, suction 
specific speed is defined as 


S = N /q/(NPSH) 3/4 


(3-1) 
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DIMENSIONLESS CAVITATION COMPLIANCE (C b /Ai) 10 



CAVITATION NUMBER RATIO, K/K* 


Fig. 8. Dimensionless Pump Cavitation 
Compliance Versus K/K* 
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which is a dimensionless parameter but, by convention, in engineering units 
Q is in gpm, N is in rpm and NPSH is in ft. For fixed N and Q, we have 


_S_ 

S* 


/npsh* \ 3;/4 

l NPSH J 


(3-2) 


Thus, the NPSH ratio would be an equivalent of the S ratio for correlating the 
data. In terms of K* and 0^, we can also write the S ratio as 


S 

S* 


K* (1 + 0^) + 0^ 
K (1 + 0*) + 0^ 


3/4 


(3-3) 


Figure 9 shows C^/A^ versus S/S* for all the pumps. Note that the 

J-2 fuel pump correlates better with the other Saturn pumps as compared to 

Fig. 8, although the J-2 oxidizer does not. The reason may be that for the 

? 

J-2 fuel pump K* is quite uncertain because of the dominance of 0^ over K* p 
whereas S/S* is less subject to this uncertainty [see Eq. (3-3)]. The J-2 
oxidizer pump has a unique feature in that there is about a 15 -percent 
recirculation of flow from the pump discharge to the inducer discharge. 

This recirculation pressurizes the inducer discharge and would be expected 
to inhibit cavitation somewhat. The consequences would be a lower com- 
pliance as shown in Fig. 9. 

In still another approach, we account for the position of the pump 
operating point relative to the region of heavy cavitation, where the latter is 
based on the effect of cavitation on pump performance. This region is con- 
sidered to extend from the point of head breakdown at K* to the initial point 
where the 4* versus K curve has zero slope, denoted by the cavitation number 
Kj. In this scheme we plot the dimensionless cavitation compliance versus 
the abscissa (K-K*)/ (Kj-K*) and the result is shown in Fig. 10. Note that 
when the abscissa is less than one, pump operation is within the region of 
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Fig. 9. Dimensionless Pump Cavitation 
Compliance Versus S/S * 
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Dimensionless Pump Cavitation Compliance 
Versus (K-K*) / K T -K* 



heavy cavitation and that this is true for all the Titan pumps (87 and 91). 
However, all the Saturn pumps operate above the region where cavitation 
has a substantial effect on pump performance. While this method of plotting 
seems to be helpful in correlating the Titan pumps, it is poorer for the 
Saturn pumps than K/K* or S/S* as seen on Figs. 8 and 9, respectively. 

For purposes of comparison with the various mathematical models 
postulated for cavitation compliance in the remainder of this report, the S/S 
plotting of the data on Fig. 9 will be employed. The best of the models are 
then also compared to the compliance data on the (K-K*)/ (Kj.-K*) basis. 

As an aside, it is interesting to note that there is no pattern on 
Figs. 8 through 10 relative to the TSH effect and pump size, except perhaps 
for a TSH effect for liquid hydrogen. If in fact these two effects are auto- 
matically taken into account by the nondimensionalization, the cavitation 
compliance data for a nonhydrogen pump can be obtained in advance of full- 
scale tests by testing a subscale model with water. Then the cavitation 
compliance of the full-scale pump with propellant can be predicted using a 
nondimens ional form such as that used for Fig. 9 , and cavitation perfor- 
mance data for the full-scale pump. The experimental determination of 
performance of a subscale model using water is a commonly employed tool 
in pump development, except that propellant is used for a liquid hydrogen 
pump. The dynamic testing would be merely another form of testing of the 
same development hardware. 
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4 . MATHEMATICAL MODELS OF PUMP 
CAVITATION COMPLIANCE 


The mathematical models developed for this study are based on 
empirical prediction of the total volume of cavitation bubbles as a function 
of the pump operating point. This is followed by a determination of com- 
pliance as the product of propellant density and the derivative of the volume 
with respect to pressure. We will consider three sources of cavitation: 

(1) cavitation attached to the blades assuming a uniform velocity at the 
inducer inlet, (2) additional cavitation resulting from the distortion of the 
velocity distribution at the inducer inlet owing to backflow, and (3) cavitation 
resulting from the presence of entrained gas in the propellant. 

4. 1 BUBBLE VOLUME 

The general approach to calculating the bubble volume V^ is to 
establish the bubble cross-sectional area at any radius A^(r) and then inte- 
grate with respect to radius 


r^ 

V b = f A b (r)dr (4-1) 

r h 

The free -streamline -wake theory and a number of empirical approaches are 
considered for calculation of the bubble area and volume. 

4. 1. 1 Bubble Volume from Free -Streamline -Wake Theory 

This theory is based on the work by Striping and Acosta (Refs. 5 
and 6). They employ a free -streamline -wake theory for a cascade of hydro- 
foils to calculate the cavity length c and maximum cavity height h for partial 
cavitation as depicted in Fig. 11, Reference 5 contains the resulting values 
of' h/d and c /d. 
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Fig. 11. Sketch of Partly CavLtating Cascade 


Appendix A describes the equations and method of approach for cal- 
culating the bubble volume and cavitation compliance using the free - 
streamline -wake theory. The following equation is developed for bubble 
volume as a function of cavitation number 


r =r . 




8 2 

3 




r 2 ) 


Ar 


r=r. 


(4-2) 


where c, h, and d are calculated for each radius and c is a function of 
cavitation number. 

4.1.2 Empirical Approaches to Bubble Volume 

The bubble area at any radius is assumed to be related to the bubble 
area at cavitation breakdown. Figure 12a shows a typical inducer with four 
blades. In the channel formed by blades A and B, we depict at radius r an 
area occupied by bubbles at head breakdown A£ (r)/z . The inlet span (d) is 
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filled with bubbles and the bubbles fill the channel (shown on the suction side 
of the blade A between c^ and c^ and on the pressure side of blade B between 
c ^ to c^. The shape of the area is quite complex. For simplicity, the actual 
bubble area between two blades is represented by an equivalent area as 
shown in Fig. 12b, where d is the span and c is an effective chord length on 
the suction side of the blade and a parabolic curve forms the third side. 
Mathematically, the area of bubbles at breakdown, as a function of the radius 
r, assuming c = d (solidity of one), is then 


A b (r) 


2 (2Trr) 2 

3 z 


(4-3) 


An alternative assumption is that the chord length of the bubble c^ at break- 
down is the product of the stay time of the bubble and the relative velocity 

c = t W(r) 

where c = chord length of bubble, t = stay time, and W(r) = relative velocity 
at pump inlet at any radius. Then 

A*(r) = | (2-rrr ) t W(r) (4-4) 

The stay time t will be determined empirically and W(r) is a function of 
radius. Next, it is assumed that the area of bubbles at radius r, for a 
cavitation number above breakdown, is equal to the product of the area of 
the bubble at breakdown times any of the ratios K * / K , S/S*, or (Kj-K*)/ 

( K - K * ) raised to the exponent n 

V r > = A b ,r| (irf < 4 - 5 ' 

or 

A b ( r ) = A^(r)(gj) (4-6) 
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or 


n 


A b (r) = Aj(r) 


Kj - K* 
K - K* 


(4-7) 


Either Eq. (4-3) or Eq. (4-4) can be employed for A£(r). 

4. 2 BLADE CAVITATION COMPLIANCE - 

MATHEMATICAL MODEL 

Our assumption for determining the compliance due to cavitation 
attached to the blades is that the velocity across the inlet is uniform. All the 
methods of determining bubble volume described in Section 4. 1 are applicable. 
The free -streamline -wake theory is a theoretical one, based on potential 
flow theory. This theory can be modified empirically by accounting for the 
difference between the theoretical cavitation number at breakdown (Ref. 5) 

2 sin a. sin ((3. -a. ) 

K * (theory) - p— - — (4-8) 

' y f 1 I cos p. 

r L 

and the value determined by test, as given in Table I. We can then define a 
modified K for use with the free - stre amline -wake theory, called 
such that 


^ mod __ K* (theory) _ j, 

K K* (test) ratio 


The use of K „ tends to correct for the effects of losses, blade thickness, 
mod 

and TSH effects. 

4. 3 BACKFLOW CAVITATION COMPLIANCE THEORY - 

MATHEMATICAL MODEL 

4. 3. 1 General Discussion of Backflow 

Backflow in an impeller is defined as the secondary flow that travels 
in an upstream direction at the impeller inlet. In a previous study (Ref. 2), 
an empirical method was developed for predicting cavitation compliance for 
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the unshrouded impellers employed in Titan pumps. In that approach it was 
assumed that the impeller tip clearance of an unshrouded impeller (radial 
clearance between the impeller tip and the housing from the inlet eye to the 
discharge of the impeller) provided a path for backflow and that the degree 
of backflow is affected by impeller performance. The goal of Ref. 2 was to 
investigate empirically the importance of the backflow in a cavitation regime 
using the tip clearance as a path. For the shrouded impeller, however, a 
different approach is needed. In the following discussion, a theory is pro- 
posed for backflow in shrouded and unshrouded impellers. 

Backflow takes place as a result of the normal blade -to -blade rela- 
tive circulation within the impeller channels, as explained in Ref. 7 and 
depicted in Fig. 13. The fluid moves relative to the blade from the suction 
side to the pressure side around the discharge side of the impeller (opposite 
to the direction of rotation). It then travels on the pressure side of the blade 
towards the inlet (opposite to the direction of through flow), and then turns 
from the pressure side to the suction side at the inlet (in the direction of 
rotation). Such a relative circulation can separate at the inlet and travel 
upstream of the pump. The relative circulation is maintained entirely within 
the channel if the rate of pressure increase along the blade is small (that is, 
if the ,! blade loading" is low). With a high blade loading it becomes possible 
for the high velocity, high momentum flow traveling along the pressure side 
to separate before turning to the suction side and then travel upstream of 
the impeller. In our view, the backflow is an extension of the normal blade - 
to-blade relative circulation required for the energy transfer from the blade 
to the fluid. With this concept, backflow could occur in either shrouded or 
unshrouded impellers at either the eye or hub of the inducer, depending on 
the condition of blade loading. In general, a design having a high flow co- 
efficient, a high specific speed, a low blade-angle turning rate, and a large 
number of blades will tend to have low blade loading and, as a result, the 
backflow' will be low. 
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In addition to backflow caused by channel recirculation, an un- 
shrouded impeller can exhibit a secondary crossflow traveling across the 
blade tip from the pressure side to the suction side. Depending on the blade 
loading and on the ratio of tip clearance to the impeller channel span, this 
crossflow could initiate or amplify a backflow within the impeller channels. 

In the shrouded impeller, the channel is closed and such crossflow cannot 
occur . 

Backflow and the interaction between backflow and cavitation is 
evident in the test data of the LR87 fuel pump as presented in Ref. 8. The 
pump tests were conducted with a series of static pressure taps located cir- 
cumferentially about the pump inlet and extending about one pump-inlet dia- 
meter upstream of the pump. The test data indicated that, at the low suction 
specific speed values, the pump-inlet pressure was several times higher 
than the normal pump-inlet pressure, indicating a strong backflow. An 
example is shown in Fig. 14. As the suction specific speed was increased, 
the pressure rise was reduced and finally at cavitation levels close to head 
breakdown the pressure rise at the pump inlet disappeared. This indicates 
an interaction between cavitation and backflow. Let us assume that the back- 
flow energy traveling upstream with high momentum is independent of the 
degree of cavitation. Then at low cavitation levels this energy is indicated 
by a high pressure rise at the pump inlet. As the suction pressure is reduced 
and cavitation is increased, the backflow momentum retards the bubble 
forming process, delaying cavitation, and as a result loses part of its 
strength. Finally, at deep cavitation the backflow momentum is absorbed 
by the strong cavitation and the inlet pressure is constrained by the vapor 
pressure. Therefore, a high backflow could lead to a reduction or delay 
in pump cavitation. This suggests that the strength of the backflow can be 
measured by the magnitude of the cavitation number at head breakdown. 

In support of the foregoing explanation. Ref. 9 shows a comparison 
between cavitation performance of a shrouded and an unshrouded impeller. 

The pump is unique, having a high flow coefficient, a low specific speed. 
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Fig. 13. Relative Circulation Within an Impeller Channel 


PUMP SPEED - 133 rev/sec, 8000 rpm 
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Fig. 14. Indicated Backflow in the LR87 Fuel Pump 
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PRESSURE ABOVE UPSTEAM VALUE 



and a low ratio of tip clearance to the impeller channel span- -all indicating 
a poor blade loading. This type of pump will have a high backflow whether 
it is shrouded or not. The test data indicated that the cavitation performance 
improved (i.e., lower cavitation number at head breakdown) with the un- 
shrouded impeller and improved even more when tip clearance was increased. 
A similar performance is also reported in Ref. 10 testing several pumps 
in a wide range of specific speed. These results correlate with the fore- 
going postulated inhibiting effect of backflow on cavitation. 

4. 3. 2 Pump Performance Related to Backflow 

As was discussed in the preceding section, backflow is a function 
of pump performance. To indicate the relative performance of the pumps 
under investigation, an attempt is made to calculate performance and design 
parameters important to backflow. Some of these parameters such as 
specific speed, efficiency, flow coefficient, cavitation number at breakdown, 
head coefficient and total head coefficient are already tabulated in Table I. 

The calculation method for other parameters such as impeller solidity and 
blade loading are now discussed and the calculated values are summarized in 
Table II. 


The method described in Ref. 5 is used to calculate the solidity t 
of an impeller, yielding 


T 


Z 


In (r /r ) 
tm lm 

2 sinp 


(4-10) 


where the equivalent number of blades z is derived from 

eq 

z 

Z “ — — 

eq smV 


(4-11) 


As shown in Table II, all Titan pumps and the H-l pumps have solidities of 
over 2. The new Rocketdyne pumps, F-l and J-2, have relatively low soli- 
dities. 
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Table II. Pump Performance Prediction Parameters 
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The diffusion factor A is a measure of blade loading in the impeller. 
As given in Ref. 11 


W t 0. 75 + Th 



/ U, \ 

U. 1 

< w l /0 t> 

£ U - — 

r\ u ./ 

) + X 


(4-12) 


where 4 1 is the actual Euler head coefficient given by 


4 Th = + oo 


(4-13) 


The ideal Euler head coefficient 4* is 

oo 




u t c v¥"ir t 


oo 


(4-14) 


and p is the slip factor. The slip factor expresses the effect of the deviation 
of relative flow velocity at the trailing part of the suction side of the blade 
caused by the relative circulation discussed in Section 4.3. 1. An equation 
developed in Ref. 13 is used for the slip factor 


P = 1 - 


2. 4 sin Y 
z -I 3 


1 + 


0. 833 


(4-15; 


A high value indicates a highly loaded impeller. As shown, all Titan pump 
blades are highly loaded. Among the Rocketdyne pumps the J-2 oxidizer and 
F-l fuel are relatively light in loading, whereas the H-l pumps and F-l 
oxidizer pumps are high in loading. 

4.3.3 Pump Performance Prediction 

Table II contains the actual Euler head coefficient for the pumps. 
This head coefficient expresses the torque transmitted from the blade to the 
fluid. In addition to the Euler head coefficient, there are external losses 
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such as disc friction, bearing losses, and secondary flow losses which, when 
added to the Euler head, result in the total head coefficient of the pump, The 
total head coefficient 4^, expressing the torque input to the shaft, is related 
to the head coefficient ^ and the efficiency rj by 




T 



(4-16) 


But, the total head coefficient is the sum of the Euler head coefficient and a 
term accounting for the external losses 


4 ^ 


T 



+ 


Ext 


(4-17) 


Note that if 4^ is calculated correctly, 4* always has a greater value than 
4^. An important parameter in 4^^ ls U^Cu^/g, where Cu^ is the induced 
prewhirl at the pump inlet. For the values tabulated in Table II, it is as- 
sumed that the prewhirl at the pump inlet is zero. There is a possibility that, 
with a strong backflow, a prewhirl is induced in the direction of rotation at 
the pump inlet resulting in a lower value of 4 J f p^ than tabulated in Table II. 

It is also important to note that the ratio of U^/U^ ^ or * n< ^ ucer e Y e iuipeller 
discharge diameter ratio) is an important parameter in determining the sig- 
nificance of the prewhirl component of the Euler head [see Eq. (4-14)]. The 
higher U^/U , the more effective is the prewhirl contribution. Therefore, 
the values of U^/U for all the pumps are tabulated in Table II. 

4.3.4 Backflow Flow Rate 

The backflow flow rate is a function of the impeller performance 

13 

parameters that were just discussed. The most important parameters 
influencing backflow are the impeller diffusion factor A, the total head co- 
efficient and the specific speed Ng. The ratio L of backflow flow rate 

to the incoming flow rate Qq is 

G> A^rp 

L = B = (constant) _I (4 .1 8) 

u o s 
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The value of the constant is based on the assumption that a good pump will 
have the following characteristics: A = 0. 6, 4^ = 0. 56 = 0.45 with 

r\ ~ 0. 80), N q = 0. 093 (1600 in engineering units), and L = 0. 0Z. The result 
is 


L 


5. 55 x 10 



(4-19) 


4. 3. 5 Velocity Distribution with Backflow 

To determine the velocity distribution at the inducer inlet in the 
presence of backflow, it is assumed that flow at the tip of the inducer travels 
upstream. The effect is to reduce the through-flow area. Based on the work 
presented in Refs. 8 and 1Z, Fig. 15 indicates the character of the velocity 
distribution. Three regions of flow are indicated: backflow, transition and a 
normal constant velocity region. It is assumed that the velocity profiles for 
both the backflow and transition regions are parabolic and that the radial 
height of both regions is the same. A detailed analysis of the velocity distri- 
bution and other related backflow parameters is given in Appendix B. 

4.3.6 Cavitation Compliance with Backflow 

The velocity distribution defined in the previous section is used to 
establish the relative velocity at each radius. For example, see Appendix B 
for a diagram of the velocity triangle at the tip of the inducer eye. Given the 
relative velocity, the cavitation number can be determined at each radius. 
Then any of the three empirical methods for establishing the bubble area 
A^(r), as given by Eqs. (4-5) through (4-7), can be applied. Then the total 
bubble volume is found by integration of the bubble area with respect to 
radius. This is accomplished for a range of cavitation numbers. Finally, 
the cavitation compliance is found by differentiating the bubble volume with 
respect to pressure and multiplying by the fluid density. 
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Assumed Backflow Velocity Distribution at the 
Inducer Inlet 


4.4 


GASEOUS CAVITATION - MATHEMATICAL MODEL 


Dissolved or entrained gases in the propellants tend to increase the 
bubble volume and compliance. We shall only consider entrained gases. In 
Appendix C, a model is developed for the volume of entrained gas within the 
cavitation region, V 

G 


G 


4 2.3 3 

= 3 " (r. - r ) 


RT P 
S k L 


X 


G 


( 4 - 20 ) 


The total compliance is then the result of the sum of the cavitation and 
gaseous volumes 


C 


b 


d < v b + V 

dP 

s 


( 4 - 21 ) 
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5. EVALUATION OF MATHEMATICAL MODELS 


Results of the mathematical models for blade cavitation and back- 
flow were compared with the experimental results. Results are first pre- 
sented for the assumption of a solidity of one for the bubble chord length at 
breakdown [Eq. (4-3)]. Results for bubble chord length based on stay time 
[Eq. (4-4)] are then presented for a selected number of the mathematical 
models. For the gaseous cavitation model no experimental results were 
available; a sample case for each of the Saturn pumps is presented to show 
the influence of entrained gas content in the propellant on the compliance. 

5. 1 FREE-STREAMLINE-WAKE THEORY FOR BLADE 

CAVITATION COMPLIANCE 

As stated in Section 4.Z, the Stripling- Acosta free -streamline -wake 

theory can be used for the mathematical model of blade cavitation compliance. 

The J-Z fuel pump, being an axial pump, is excluded from the investigation. 

Figures 16 and 17 present the cavitation compliance for the F-l, H- 1 and 

J-Z pumps (except for the J-Z fuel pump) showing the experimental results 

and also the predicted values. Figure 16 shows the predicted values without 

modification of the cavitation number given by Eq. (4-9) (that is, ^ rdi ^ lo = *)• 

Figure 17 employs a modified cavitation number, with the values of ^ Vdi ^ 0 

presented in Table III. An examination of values of K appearing in 

Table III shows that for all the pumps, because of the losses and real fluid 

effects, the values of K are less than one. In addition, because of the 

TSH effects of the oxidizer, K ,. values for the oxidizer pumps are higher 

ratio 

than K ra ^ 0 for the fuel pumps. Therefore, the use of seems to pro- 

vide a logical basis for blade cavitation compliance, since losses and pro- 
pellant effects are automatically taken into account. 

Figure 16 shows that the prediction method using the free -streamline - 
wake theory yields much lower values than the test results. In the extreme 
case, the predicted H-l fuel compliance is two orders of magnitude lower. 
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Fig. 16. Test Versus Blade Cavitation Prediction Using 
Free -Streamline -Wake Theory 
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TEST DATA 

BLADE CAVITATION 

ASSUMPTIONS: 

1. VAPOR VOLUME BASED ON 
STRIPLING AND ACOSTA 
FREE-STREAMLINE-WAKE 
THEORY 



SUCTION SPECIFIC SPEED RATIO, S/S* 


Fig. 17. Test Versus Blade Cavitation Prediction Using 
Free -Streamline -Wake Theory With ^ 
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With use of the values, the predicted compliances shown in Fig, 17 

increase for some of the pumps, but the values are still quite low. We con- 
clude that the free -streamline -wake theory does not apply to the cavitating 
regimes of rocket pumps. It is probably limited to non- cavitating or incipient 
cavitation region. 


Table III. Values of K 

ratio 


Engine 

F- 


J-Z 

H- 

I 

Pump 

OxLdizer 

Fuel 

Oxidizer 

Oxidize r 

Fuel 

K 

ratio 

0. 687 

0. 307 

0. 795 

0. 785 

0. 132 


5. Z RESULTS FOR CAVITATION COMPLIANCE 

CONSIDERING BACKFLOW 

By the method presented in Section 4.3.6, the cavitation compliance 
considering the effect of backflow on the inlet velocity distribution is calcu- 
lated for all the Rocketdyne pumps except the J-Z fuel pump. The unity 
solidity assumption [Eq. (4-3)] is used for these results. In an initial trial, 
a cavitation number ratio [Eq. (4-5)] with the exponent n = 1/Z is used for 
the backflow region, and the free -streamline -wake theory is used for the 
transition and normal regions. The computed results are shown in Fig. 18. 
As can be observed the predicted values are low. 

For a next trial, Eq. (4-5) with n = 1/Z is used for all three flow 
regions (backflow, transition and normal). At the same time, to emphasize 
the effect of incidence angle, a modified value of relative velocity W^ is used 
in calculating the cavitation number. It is assumed that, due to free- 
stream diffusion, the relative velocity at any radius entering the blade is 
determined as follows: 


W, (r ) = C /sin p 
1 m r 


(5-1) 
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TEST DATA 

BLADE CAVITATION 


ASSUMPTIONS: 


1 


THE TRANSITION AND NORMAL 
SECTIONS VAPOR VOLUME 


CALCULATION BASE ON THE 
STRIPLING AND ACOSTA FREE- 
STREAMLINE-WAKE THEORY 


WITH K 


MOD 



2. BACKFLOW BUBBLE VOLUME 
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Fig. 18. Test Versus Backflow Cavitation Prediction With 
(K*/K)°- 5 and r = 1 
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where p equals the blade angle when the incidence angle a <4 deg. For 
values of a > 4 deg, p equals the approach relative velocity angle (3^ plus 
4 deg. It is further assumed that the shock loss at the inlet is negligible. 
Recall that Wj is used in the calculation of K [Eq. (2-1)]. These changes 
were incorporated in the computer program and the results presented in 
Fig. 19- As shown, the predictions are now more comparable with the data, 
although the calculated compliances have higher slopes than do the test 
results. The calculated compliances for the F-l pumps cross the test 
results, the H-l pumps have low calculated compliances, and the J-2 oxi- 
dizer pump has high calculated values. As discussed in Section 2, the J-2 
oxidizer pump is expected to have a higher predicted value because of the 
effect of recirculation to the inducer exit in inhibiting cavitation. 

For a third trial the suction specific speed ratio [Eq. (4-6)] is used 
with n = 1/2. In terms of an NPSH ratio, related to the S ratio by Eq. (3-2), 
the bubble area at any radius is 

A b (r) = I V r2 ( npsh */ npsh ) 3n/4 (5 - 2) 

The results of the compliance calculations are shown in Fig. 20. Comparing 
Fig. 20 with Fig. 19, some improvement is seen by using the suction 
specific speed ratio. 

To observe the effect of a change of the exponent, n was varied 
from 1/2 to 0. 15. The best match was obtained at 0.2 and these results 
are shown in Fig. 2 1. The compliances for the J-2 oxidizer pump, the F-l 
pumps, and the H-l oxidizer pump are closer to the test results, but the 
compliance for H-l fuel pump is slightly lower than the data presented in 
Fig. 20. 
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Fig. 19. Test Versus Backflow Cavitation Prediction 
With (K*/K)°- 5 and r = 1 
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Fig. 20. Test Versus Backflow Cavitation Prediction 
With (S/S*)0- 5 and t = 1 
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DIMENSIONLESS CAVITATION COMPLIANCE |C b /A]| x 10 


TEST DATA 


BACKFLOW CAVITATION 



Fig. 21. Test Versus Backflow Cavitation Prediction 
With (S/S*)°* 2 and t = 1 
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RESULTS OF EMPIRICAL APPROACHES TO BLADE 
CAVITATION COMPLIANCE 


5. 3 


In the previous section, we accounted for the effect of backflow on 
the radial distribution of inlet velocity. We considered both the cavitation 
number ratio and the suction specific speed ratio as a basis for the bubble 
area at each radius. We took the cavitation number to vary with radius 
because it depends on relative velocity. However, in our treatment, the 
suction specific speed ratio is velocity independent and thus is uninfluenced by 
backflow. For the following reasons we are led to the position that the 
consideration of backflow is not warranted. First, the use of suction specific 
speed ratio which is independent of backflow led to as good a correlation with 
the test data as did the cavitation number ratio, which is dependent on back- 
flow. Second, the consideration of backflow is a significant complication in 
computations and requires a number of assumptions. Therefore, we now 
restrict ourselves to the computation of blade cavitation compliance, which 
implies a uniform velocity at the inducer inlet. 

Using S/S* and unity solidity based [Eqs. (4-3) and (4-6)], we find 
the bubble volume to be 


V, 


8 ir 2 / 3 3\ /NPSH*\ 3n/ ^ 4 

9 z \ r i " r h/ V NPSH / 


(5-3) 


By definition NPSH is related to the static pressure P g 
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Thus Eq. (5-3) becomes 
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Taking the derivative with respect to P 
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In dimensionless form 
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Using Eq. (5 - 5) 
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In engineering units 
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In terms of cavitation number 
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In engineering units, the only change in Eq. (5-11) is a change of the coeffi- 
cient from 4/3 to 1/9. 

Figure 22 shows the results of calculations for all the pumps 
(including the J-2 fuel and the Titan pumps) using Eq. (5-9) with n = 0. 2. 
Note that the results in Fig. 21 are identical to the corresponding results on 
Fig. 22. The prediction for the J-2 fuel pump is low compared to the test 
results. The Titan results except for LR91 fuel show the correct relative 
position, but have a flatter shape. 

If, instead of the unity solidity assumption, we employ the stay 
time basis for bubble area [Eq. (4-4)], the following equations result: 
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By using Eq. (5-5) to substitute for NPSH 
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TEST DATA 

BLADE CAVITATION 



0.1 0.2 0.3 0.4 0.6 0.8 1.0 


SUCTION SPECIFIC SPEED RATIO, S/S* 

Fig. 22. Test Versus Blade Cavitation Prediction With 
(S/S#)°* 2 and r = 1 


- 51 - 



and 
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In dimensionless form 
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In terms of NPSH 
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In engineering units the coefficient ^g 11 is omitted from Eq. (5-16). Finally, 
in terms of cavitation number 
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(5-17) 


Equation (5-17) is used in the computer program where W(r) is determined 
for each radius by the method described in Section 5.2 using Eq. (5-1) with 
free stream diffusion at the inlet of the inducer. It is interesting to note that 
number of blades, z, does not appear in Eqs. (5-12) through (5-17). Figure 23 
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DIMENSIONLESS CAVITATION COMPLIANCE (C b /A 1 ) 10 










shows the results of calculation using Eq. ( 5-17) with n = 0. 2 and a value 
of t = 4. 8 msec based on an empirical fit of the F-l oxidizer data. Figure 23 
shows reasonable correlation between the test data and prediction for most 
of the Rocketdyne pumps. For the Titan pumps, although the values for 
LR87 fuel are close, in general the prediction is poor. The predictions for 
F-l pumps are very close and the values for H- 1 and J-2 fuel are reasonable. 
The J-2 oxidizer pump is not predicted well. As pointed out earlier, the 
recirculatory flow may be inhibiting the cavitation of the J-2 oxidizer pump. 
All in all, Eq. (5-17) with a single empirically determined value of t appears 
to yield the best results. Perhaps a further improvement could be realized 
by accounting for the effect of the heat transfer properties of the propellant 
on the stay time t, although the propellant heat-transfer difference is not 
seen between RP-1 and LOX for the F-l pumps. 

Up to this point, all the evaluations were based on presenting the 
dimensionless cavitation compliance versus S/S*. We now use Eq. (4-7) 
to obtain the compliance as a function of (Kj - K*)/ (K - K * ) . Both the unity 
solidity and constant stay time assumptions are used, Eqs. (4-3) and (4-4), 
respectively. 

The results are shown in Fig. 24 for the unity solidity assumption 
and in Fig. 25 for the constant stay time assumption. In both cases, n = 0. 2 
is assumed. Figure 24 shows that except for H-l fuel and J-2 fuel the pre- 
diction for all Rocketdyne pumps is reasonable. The prediction for the Titan 
pumps is poor. Figure 25 exhibits the same trend as shown in Fig. 23 for 
the Rocketdyne pumps, but the values are somewhat closer to the test results. 
Again, the H-l fuel and J-2 oxidizer are predicted at low values. However, 
at low abscissa values the H-l fuel prediction is close to the test results. 

For the Titan pumps it is interesting to note that the prediction method shows 
higher values but the predicted values cross the same way that the test 
results do. It is postulated that the bubble stay time of 4. 8 x 10 sec 
(based on the F-l oxidizer pump) does not apply to the Titan pumps. The 
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Fig. 24. Test Versus Blade Cavitation Prediction With 
[(K T -K*)/(K-K*)1°‘ 15 and r = 1 
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Fig. 25. Test Versus Blade Cavitation Prediction With 
[(K t -K*)/(K-K*)]°- 15 and t = 4. 8 msec 




Titan inducers are cambered with shorter blade length and the bubble length 
is expected to be shorter due to the early energy input to the fluid by the 
blade. For example, if the value of stay time is reduced by 4 times, to 
1. 2 msec for the Titan pumps, the prediction method will match the test 
results for LR87 oxidizer and LR87 fuel pumps. The same logic could be 
used for the J-2 oxidizer pump where the recir culatory flow to the inducer 
discharge pressurizes the inducer trailing edge and reduces the bubble chord 
length or stay time. For example, a stay time of 2.4 msec for the J-2 
oxidizer pump will match the prediction values with the test results. 

5.4 GASEOUS CAVITATION COMPLIANCE 

In Section 4. 4, an approach was presented for calculating the gaseous 
cavitation compliance. To demonstrate gaseous cavitation compliance, a 
sample calculation was made assuming helium gas entrained with a mass 
fraction X^, of 0. 01 percent. The results are shown in Fig. 26 for all the 
Rocketdyne pumps except the J-2 fuel. As can be observed, the gas content 
has a significant effect on the cavitation compliance. If one assumes a static 
pressure of 46 psia and a propellant temperature of 70°F for the F-l fuel 
pump, the above gas content corresponds to a 15.5 percent by volume of 
helium in RP-1. 

5. 5 CURVE FITTING EQUATION 

Although the prediction method provided a reasonable way to calcu- 
late the cavitation compliance, it could not be used successfully to predict 
the compliance for all the pumps. In an attempt to use the data directly, the 
following curve fit equation was developed to match the Rocketdyne data 

2 3 

(C^/A^ x 10 5 = - 22. 4 + 286 (^) - 274 (^) + 92. 6(Jr) (5-18) 

This function is plotted in Fig. 27. In this match of the data, both the J-2 
fuel and the J-2 oxidizer pumps are de -emphasized. 

Due to poor correlation of data for the Titan pumps, no attempt was 
made to curve fit the compliance for these pumps. 
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0.2 


0.4 0.6 0.8 1.0 


0 . 1 

SUCTION SPECIFIC SPEED RATIO, S/S* 


Fig. 26. Predicted Total Compliance With a Gas 
Content of = 0,01 Percent 












6. CONCLUSIONS 


The following conclusions have resulted from the study of empirical 
correlation of the cavitation compliance test data and from the studies of 
mathematical models for predicting cavitation compliance: 

(a) The cavitation compliance test data of the H- 1 and F - 1 Saturn 

pumps can be correlated quite well by a particular nondimens ionalization, as 
shown in Figs. 8 and 9. This indicates that the influence of pump design, 
pump size, and propellant properties are inherently taken into account. The 
nondimens ional cavitation compliance is expressed as the volumetric com- 
pliance (volume change per unit pressure change) times the propellant weight 
density, divided by the area of the pump inlet. The nondimensional per- 
formance parameter can be any of the following divided by its value at head 
breakdown: cavitation number, suction specific speed, or net positive suc- 

tion head. An expression for a best-fit curve is given by Eq. (5-11) and 
shown in Fig. 27 relative to the test data. The J-2 oxidizer pump is less 
well correlated because of a unique feature, recirculation to the inducer 
discharge, which tends to reduce cavitation. The comparison of the J-2 
fuel pump with the others is made uncertain because of a lack of accuracy in 
establishing head breakdown from the available pump performance data. 

(b) The above conclusion further suggests that, for the Rocketdyne 
family of pump designs, the method of empirical correlation can also be 
employed to relate the cavitation compliance of a subscale pump with water 

to that of a full-scale pump with propellant. A possible exception is that 
liquid hydrogen may not have a suitable replacement for this purpose. 

(c) The cavitation compliance test data of the LR87 and 91 Titan 
pumps can be correlated by the nondimensionalization shown in Fig. 10. 

The nondimensional compliance is the same as that identified in the first con- 
clusion, but the performance parameter is a cavitation-number difference 
ratio which expresses the relative position of the cavitation number within the 
range from head breakdown to the point of zero slope of the head rise versus 
NPSH characteristic. 
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(d) The blade cavitation compliance predicted by the free- 
s treamline -wake theory is too low, indicating that the assumptions of this 
theory (potential flow, zero blade thickness) render it inapplicable to the 
cavitation regiine of rocket pumps. 

(e ) Consideration of the effect of backflow on the velocity distri- 
bution across the inducer inlet provides no improvement in the agreement 
of the mathematical models with test data. 

(f) The cavitation compliance of the Saturn H-l and F-l pumps 
can be matched by assuming a bubble proportional to the 0. 2 power of the 
suction specific speed ratio times the bubble volume at head breakdown 
(see Fig. 23). An equally good alternative is the assumption that the bubble 
volume is proportional to the 0. 15 power of the cavitation-number difference 
ratio times the bubble volume at head breakdown (see Fig. 25). The volume 
at breakdown is found by extending the pump inlet area along the inducer 
blades for a chord length equal to the product of the relative velocity and a 
stay time empirically determined to be 4. 8 msec (see Fig. 23). The use of 
a chord length corresponding to a solidity of one yields a somewhat poorer 
correlation; see Fig. 22. On this basis the J-2 oxidizer pump appears to 
have a shorter stay time, perhaps the result of the recirculation flow's 
confinement of the bubbles to a shorter distance along the inducer blades. 

(g) The cavitation compliance of the Titan pumps can be matched 
by assuming a bubble volume proportional to the 0. 15 power of the cavitation- 
number difference ratio times the bubble volume at head breakdown; see 

Fig. 24. The volume at breakdown is like that in conclusion (f), except that 
the stay time is approximately 1.2 msec. A shorter stay time relative to 
the Rocketdyne pumps is to be expected because of the cambered design of 
the inducer blades of the Titan pumps. 

(h) A method has been formulated for predicting the contribution 
of entrained gas to cavitation compliance. No data are available for evalua- 
tion of the method. 
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7. RECOMMENDATIONS 


The following recommendations are made for prediction of cavitation 
compliance for new pump designs in the Rocketdyne family: 

(a) Prior to the availability of experimental data on a subscale or 
fulls cale pump, it is recommended that two methods be employed to estimate 
the nondimensional cavitation compliance. They are: 

1. The best-fit curve of Eq. (5-18) shown in Fig. 27, and 

2. The mathematical model which assumes a bubble volume 
proportional to the 0.2 power of the suction specific 
speed ratio times the bubble volume at head breakdown, 
along with a stay time of 4. 8 msec, as exemplified in 
Fig. 23. 

Any difference between the two estimates is an indication of their uncertainty. 
The cavitation compliance can then be estimated using the predicted position 
of head breakdown. 

(b) If subscale pumps are evaluated during the development of 
the new design, pulsing tests should be performed to establish a nondimen- 
sional plot of cavitation compliance such as that of Figs. 8 or 9. Particular 
attention should be paid to establishing the position of head breakdown during 
the performance testing. The cavitation compliance of the fulls cale pump is 
predicted by assuming it possesses the identical nondimensional cavitation 
compliance characteristic, even if water has been substituted for the pro- 
pellant for the subs cale tests. 

(c) It is recommended that the cavitation compliance of subscale 
models of the existing Rocketdyne pumps be determined experimentally to 
obtain an additional basis for empirical correlation. 
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APPENDIX A 


BLADE CAVITATION COMPLIANCE 


In Ref. 5, the free - streamline -wake theory is used to calculate the 
cavity height and length at an inducer leading edge radius 


h/d = sinp 
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where 


R = w i / w ^ 

a = P- P F 

d = blade spacing 


From Ref. 4, is given by the following equation 
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where 


w = w. \/l + K (A-4) 

c 1 

In using these equations it is assumed that 

(3 = f(r) 

In most cases, the leading edge is a radial element or has a constant 
pitch. For the radial element blades 

r tan (3 = Constant (A- 5) 

Another important parameter is the axial velocity (C m ^) flow coeffi- 
cient (0^) distribution at the leading edge as a function of radius. 

C 

m 

0 i = -IT 1 (A_6) 

i 

It is assumed that 0^ is constant along the leading edge of the reducer. 

For the initial estimation it was also assumed that the blade thickness 
was zero, and axial velocity and pressure distribution at the inlet were uni- 
form. With these assumptions in mind, the inducer blade radii from eye to 
the hub were divided into several stations and the cavity height and lengths 
were determined for each station. By using the cross-sectional area of the 
cavity and the radius increment (Ar) between the radial stations, the bubble 
volume was then determined. The following equations were used. 
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Main inputs are 


K. , r . , U. , and 0 , 

ill 1 


By definition the local cavitation number at any radius is 
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NPSH and C are constant at the inlet but W. varies with radius and assuming 
m j 

no prewhirl at the inlet 
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Equation (A-7) can be written in the following form 
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Therefore, K. can be described in terms of K. 

J 1 

K (1 + 0 Z ) 

but 

U. 

— 2 . = 

U . 



(A - 1 1 ) 


(A- 12 ) 


Therefore 
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With the assumption of no prewhirl the relative angle of flow at any radius is 


tan S = C /U. 
r f . m 1 
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using Eq. (A- 1 2 ) 
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The incidence angle 


a . 
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(A- 1 6 ) 


With the above information the values of c/d and h/d at any radius can be 
determined. The area of the cavity, for each blade, at any radius, assuming 
parabolic cavity shape is given by 
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The cavity volume along the radius for all blades 
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To obtain -dV/dP s values, 
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s 


using cavitation number equation 
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The value of dV/dPis obtained graphically by plotting V versus P g . Finally, 
the dimensionless cavitation compliance C A^ = (-dV/dP) P^/A^. 
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APPENDIX B 


BACKFLOW CAVITATION COMPLIANCE 


Figure 14 shows the assumed velocity profile for a pump with backflow 
and the profile is divided into three regions: 

1. Backflow region 

2. Transition region 

3. Normal region 

In the following discussion the velocity distribution and the related 
values required for the calculation of c/d and h/d using the free-streamline - 
wake theory are calculated. 


I. BACKFLOW REGION 


It is assumed that the backflow Euler head at the eye of inducer equals 
half of the power head generated by the pump 


U B C Ve 
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where and Cu n are determined for each radius and blade angle, 
jd r> 

Assuming that backflow follows the blade, the velocity triangle for the 
backflow is shown in Fig. 28. In Fig. 28, C^ is the absolute velocity of the 


backflow, Wg is the relative velocity and Ugthe tangential velocity, 
the meridional or axial component of absolute velocity. 
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and from the velocity triangle 
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Knowing and p^ the values of 6^ and Wg can be determined. 



Fig. 28. Velocity Triangle for the Backflow Region 


It l s assumed that C m distribution is parabolic along the radius with 

the maximum C at the outer diameter and zero value at the minimum dia- 
m 

meter as shown in Fig. 15. Thus 


C 



m Bi 

(r.-r. ) 
v i la 



(B-5) 


- 69 - 



where Cm^ corresponds to the velocity at r. Knowing the value of 
to determine the value of r- the following integration is used: 

1 U 
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At this point, the region of backflow extends from r^ to r^. And at 
every r the value of C m ^, p, U, W are known* Assuming the flow follows 
the blade, a = 0 and p = p^. It is also assumed that the static pressure at 
the backflow region equals stagnation pressure at the normal region of the 
pump. Therefore 
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where P g and C m refer to static pressure and axial component of velocity 
"normal region” which will be considered later. 
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The only variable in the above equation is Wg which is a function of 
radius, and is obtained from the velocity triangle Wg = C m g/sin (3. 
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At this point the input values are known and the empirical equations 
presented in the text are used to calculate compliance. 

II. TRANSITION AND NORMAL REGIONS 

As shown in Fig. 15, these two regions constitute the through flow 
velocity profile. The transition region extends a radial length equal to the 
radial length of the backflow region. 
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and has parabolic velocity distribution. At the upper radius the value of 
velocity is zero and at the lower rac 
velocity for the normal region). So 
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In the "normal region" 
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Combining Eqs. (B-14), (B-15) and (B-16) 
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Therefore, the C distribution is known and using the method explained 
m 

in .Appendix A the values of v Trans + N are determined for a radius range 

of r. to r, . 
la h 


The total vapor volume for a pump with backflow is 
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and -dV/dP^ is determined graphically. Finally, = -dV/dP (Pj^). 
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APPENDIX C 


GASEOUS CAVITATION COMPLIANCE 


It is assumed that gas is entrained in the propellant and the propellant 
flow is a mixture of gas and liquid. The gas properties and the mass frac- 
tion of gas in liquid are known. It is also assumed that the gas volume 
change follows perfect gas laws and the process is isothermal. It is also 
assumed that X , the mass fraction of gas in liquid, is known 

x G = m q /m l (C-l) 
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Assuming that up to a solidity of one the static pressure is basically equal to 
the suction pressure but beyond that the pressure increases and the gas 
volume becomes insignificant. At any radius 

AV = (2-rrr)^ (C-4) 

L z 


- 73 - 



using Eq. (C-3) and substituting AV T for V and AV for V„ 

-L* J-j Ci Lj 


2 Z R T p T 

A V G = 4 ~~z -p LilX G (Ar) 


(C-5 ) 


Integrating with respect to radius 
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Then total cavitation compliance 
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